Orbital Correlations in the Ferromagnetic Half-Metal CVO2 
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We deduce a model relevant for the description of the ferromagnetic half-metal Chromium dioxide 
(CVO2), widely used in magnetic recording technology. The model describes the effect of dynamical, 
local orbital correlations arising from local quantum chemistry of the material. A finite temperature 
solution of the model in d — 00 provides a natural explanation of the optical response, photoe- 
mission, resistivity and the large Woods-Saxon ratio observed in experiments. Our study confirms 
the important role of many body dynamical correlation effects for a proper understanding of the 
metallic phase of CrOi. 



PACS numbers: 75.30.Mb, 74.80.-g, 71.55.Jv 

Colossal Magnetoresistive (CMR) materials have been 
the focus of renewed theoretical and experimental inves- 
tigations [0J^] in recent years, given their obvious tech- 
nological potential. Interest has mainly been centered 
around the manganites, where the interesting ferromag- 
netic metallic regime is driven by the double-exchange 
(DE) mechanism, though a unified understanding of the 
correlated nature of the metallic state itself, as well as the 
high-T insulating paramagnet, requires consideration of 
orbital correlations and the Jahn- Teller (JT) distortion 
on an equal footing with the DE ferromagnetism. 

Some attention has also focused on Cr02, widely used 
in magnetic recording tapes. In contrast to the man- 
ganites, stoichiometric CrOi is already a ferromagnetic 
metal. Given the formal 4+ valence state of Cr, the two 
3c? electrons occupy the t2 g orbitals. One would intu- 
itively expect to form S = 1 spin on each site, and an 
antifcrromagnetic Mott insulator. Why CrOi is a ferro- 
magnetic metal instead, has been answered by Korotin 
et al H , who have carried out insightful (LDA + U) cal- 
culations for this material. Their main conclusions are: 

(i) the O 2p band(s) act, atleast partially, as electron 
(or hole) reservoirs resulting in Cr being mixed- valent 
(like Mn in doped manganites), explaining metallicity, 

(ii) an almost dispersionless majority spin band of pre- 
dominantly d character at about leV below Ep over a 
large region of the Brillouin zone. This corresponds to 
strongly localized xy orbitals completely occupied by one 
majority spin electron. On the other hand, the d states 
of predominantly d yz + zx character hybridize with the O 
2p band and disperse, crossing Ep. The Hund's rule cou- 
pling between the localized d xy spin and the spin density 
of the band d yz+zx electrons polarizes the latter, giving a 
ferromagnetic state via the double-exchange mechanism. 
Thus, both the metallicity and ferromagnetism are cor- 
related well with each other. 

A closer examination reveals that the metallic state of 
CrOi is strongly correlated, implying that many-body 
correlation effects beyond the local-density approxima- 
tion (LDA) |U (or its variants) need to be considered. 



A number of experimental observations tend to support 
such a picture: 

(1) Polarization dependent XPS measurements reveal 
substantial ligand orbital polarization. An exchange 
splitting energy of A ex ^ sp i ~ 3.2 eV was deduced [|J, 
implying substantial correlation effects. On the other 
hand, LSDA calculations yield A ex _ sp i ~ 1.8eV! 

(2) The resistivity shows a behavior characteristic of 
correlated Fermi liquids §: p(T) = p Q + AT 2 + BT 7 / 2 , 
the last term coming from the carrier scattering of 
two-magnon fluctuations in a double-exchange ferromag- 
net |tJ . A similar kind of behavior has also been observed 
in the CMR manganites in the FM metallic phase JjJ . In 
fact, CVO2 falls into the "bad metal" classification, with 
high-T (T > T C FM = 390K) resistivity exceeding the 
Mott limit |pjo| . The coefficient A is large, and in fact 
the ratio A/7 2 (7 is the coeeficient of the linear term in 
the electronic specific heat) is close to the value expected 
for heavy-fermion metals (^] , implying substantial corre- 
lation effects in the Cr <i-band. 

(3) Additional evidence comes from the optical conduc- 
tivity data, which reveals a Drude conductivity at low 
energies, followed by a broad bump around 0.8eV and 
high-energy features centered around 3eV || . While the 
Drude part is understandable within the LSDA approach, 
the other features mentioned above are seen at energies 
that are systematically 10 — 20 percent lower than those 
predicted by the LSDA calculation, showing up the im- 
portance of correlation effects in this material Q . 

(4) Additionally, recent magnetotransport measure- 
ments show interesting features. Firstly, the magnetore- 
sistance is linear in fields upto 1 Tesla, and shows no 
evident relation with the magnetization j[o|. On the 
other hand, it was found that the longitudinal and trans- 
verse MR is negative and increases linearly with field for 
T > 200K, while having a concave shape as a function 
of field for lower T |10|, passing through minima as func- 
tions of field around T ~ 200K. It is interesting to notice 
that similar behavior is also observed for manganites, 
Lai- x Sr x Mn0 3 with x = 0.175,0.2 O. An under- 
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standing of the observed crossovers around T — 200K 
requires the consideration of additional mechanism(s) 
over and above the double-exchange interaction, since 
these features occur deep inside the ferromagnetic metal- 
lic phase (T c = 390K). 

(5) Existing photoemission data on CrC>2 show fea- 
tures that are more characteristic of a semiconductor 
with vanishing DOS at the Fermi level than of a good 
metal completely consistent with the "bad metal" 
classification made in the resistivity measurements. Such 
features, also observed in the metallic FM state of the 
CMR manganites, have been linked to the strong Jahn- 
Teller distortions known to exist for CMR mangan- 
ites ]l7|] . However, neither the static J-T distortion, nor 
the doping- induced static disorder is present in Cr0 2 , 
leaving one to look for alternative scenarios to under- 
stand the suppression of low-energy spectral weight in 
PES. 

We start by noticing that the emergence of a coherent, 
strongly correlated Fermi liquid scale required to under- 
stand the above features is out of reach of LDA based 
calculations (within which a large Woods-Saxon ratio 
is inexplicable, for example), as well as of calculations 
based on pure double exchange models. This is because 
these features are observed deep inside the ferromagnetic 
metallic state, and are thus related to additional scatter- 
ing mechanisms in a half-metallic situation (where one 
deals with a fully spin-polarized band). In this paper, we 
argue that the above effects can indeed be understood 
by invoking the important role of strong local, orbital 
correlations in the t 2g sector (see below). 

An understanding of features mentioned above should 
go hand-in-hand with the basic electronic structure. In 
this context, the LSDA+U calculation carried out by Ko- 
rotin et al. || already provides the necessary ingredients 
to pursue more detailed investigations. The necessary 
ingredients seem to be: (1) Strong, local Hubbard U 
and Hund's rule coupling, Jh- (2) Hybridization of the 
dyz+xz band with the p band of the O, resulting in Cr be- 
ing mixed- valent, like Mn in the manganites. As pointed 
out in ||, the fact that Cr0 2 is in the class of negative 
charge-transfer gap materials makes this d — p hybrid 
band self-doped. The well known double-exchange mech- 
anism then drives the FM metallic phase. 

An examination of Figs. (3)-(4) of Ref. Q shows that 
the the bands crossing the Fermi level have comparable 
dxz+yz and d yz - Z x character, as well as a non-dispersive 
d xy band centered leV below \i. We shall consider the 
O band to act solely as a reservoir of carriers, serving 
to make the <i-band manifold non- half-filled @. The al- 
most complete spin polarization tells us that the local 
Hund's rule coupling is strong. If this is so, we have to 
assume that the on-site Coulomb repulsion (the intra- 
orbital Hubbard U) is stronger. Given the above, we are 
then led to consider a two-orbital Hubbard-like model in 
the double exchange (DE) limit, U, Jh >> t a ,b (a,b la- 



bel orbital indices). The relevant minimal model for our 
purpose is written as, 

H = H 1 +H 2 (1) 

where 

Hl = ~ t "H4aa c 3b^ + h.c) + U ^2{n iar n iai +a^b) 

<ij>a i 

(2) 

and 

H 2 = U ab ^ n iaan lba , - J H ^ s i ' (°f + °f) (3) 

iaa' i 

In the above, we shall assume that the effects of the 
actual bandstructure have been incorporated into the ef- 
fective hopping strengths by fitting the idealized band- 
structure to the actual one; this has been employed with 
much success for the cuprates |l3[ ]. We shall henceforth 
treat the t ab as parameters of the effective model above. 
U a b is the local, inter-orbital Coulomb repulsion, Si is 
the localized spin in the d xy orbital, and that a = d xz + yz 
and b = d yz - Z x- In the DE limit, the states correspond- 
ing to a carrier aligned antiparallel to the core-spin are 
projected out, and so one can drop the spin indices in the 
above eqn, at the price of introducing a magnetization- 
dependent hopping, tfj = <° b (S). In this situation, the 
orbital index plays the role of the spin, and so one is fi- 
nally left with a generalized Hubbard model in orbital 
space. To proceed, we will further make an assumption 
which does not qualitatively affect the final conclusions 
(see below). We will assume that the hopping matrix is 
diagonal in orbital space, tf b = tijS a i, Q|, and choose 
t aa = t\ and t bb — t 2 . Finally, we introduce new fermion 
operators, Cf = (c a + Cb)/y/2, and C| = (c a — Cb)/V^- The 
Hamiltonian now becomes, 

H = - E ^(S)(4c iCT +/i.c)+C/^n iT n a (4) 

<ij>a i 

At low-T (< T™ = 390K), the hopping is a function 
of the magnetization, £ CT (S) = t a (M) = t a [l+ < SfSj > 
/2S 2 ] 1 ' 2 with M = N- 1 £\ < Sf >, so one has to deal 
with a Hubbard model in orbital space with an orbital 
pseudo-spin and magnetization dependent hopping. In 
the above eqn. (4), if = (ii + £2), and t\ = (ti — t 2 ). 
The case t\ = t 2 has been considered in |L4|], where the 
anomalous FM metallic phase of the manganites is sought 
to be understood in terms of the non-FL state of the 
effective Falicov-Kimball model in d — 00 jl5| . An exact 
solution of the FKM in this limit reveals that the non-FL 
behavior is driven by a divergence in the local excitonic 
correlation function at low energy. In contrast, there is 
evidence that the FM metallic phase in Cr0 2 as well as 
in Lao.rSro.sMnOs is a correlated Fermi liquid at low T, 
requiring us to work with the model, eqn. (4). 
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In what follows, we will explain much of the observed 
features mentioned above in terms of the d = oo solution 
of the effective Hubbard model (eqn. (4)). To understand 
where the local FL behavior comes from, we notice that 
the model, eqn. (4) with t± = ti, can be mapped onto 
the x-ray edge problem with the recoilless local scatterer 
(the |-spins) right at the Fermi surface. In this case, 
the spectral fn. of the |-spin fermion, as well as that of 
the "exciton" jl5| is divergent at low energy. The case 
t\ 7^ t2 corresponds to allowing the scatterer to recoil, 
killing the infrared divergence and restoring FL behavior 
at low T (T < En, the recoil energy of the |-spin [lq] ). In 
the d — oo solution of the model, however, FL coherence 
sets in below the lattice coherence scale, T co h (related to 

h-h). 

It is instructive to summarize what is known about the 
Hubbard- like model, eqn. (4), in the d = oo limit [ p^[ . 
Above Tcoh, the metallic state off half- filling is not a 
Fermi liquid, and the dynamics of the carriers is de- 
termined by inelastic scattering off the effectively un- 
quenched local moments (orbital moments in our ap- 
proach). Below T co h, the local moments are screened by 
the "conduction electron" spin density, giving a Fermi 
liquid response. This coherence scale can be driven quite 
low for given t\ — ti resulting in an anomalous response 
at not too low T. In particular, the suppression of quasi- 
particle features in the local DOS for T > T co h |l6| would 
result in: (a) a non-Drude optical response, (b) a pseu- 
dogap feature in photoemission, (c) absence of the T 2 
term in the resistivity and an anomalous T-dependence 
of the Hall constant jl6| , and (d) lastly, one would ex- 
pect the FL coherent response below T co h to manifest 
itself in these responses as T is lowered. In particular, 
below T co h, one expects a T 2 component in the resistivity, 
a low energy Drude response in the optical conductivity, 
a coherent feature in the photoemission lineshape at the 
Fermi surface, and a T-independent normal Hall constant 
(the anomalous Hall effect requires more work 18|). 

Given that the para-orbital metallic state of the above 
model with t% ^ t 2 is a correlated Fermi liquid, one ex- 
pects perturbations like disorder to tilt the balance in 
favor of low-energy incoherence. In particular, with mod- 
est disorder [|1S| ] , the incoherent metallic state gives rise 
to an incoherent lineshape in PES and a pseudogapped 
behavior in the optical response, even at low-T. 

With this information in mind, we are ready to dis- 
cuss our results. In our calculations, we have employed 
a gaussian unperturbed DOS, and two values of temper- 
ature: T = 0.0113 and T = O.LD (D is the half-width 
of the gaussian). We have chosen U a b = 3-Dj, where 
Z3| is the effective "f-spin bandwidth, and the bandfilling, 
n = rif + ni = 0.8. In the d = oo approximation that 
we use, to-(M) = t a [(l + M 2 (T))/2] 1 / 2 with M(T) taken 
from the spherical model limit for the Heisenberg model. 

In d — oo, all relevant information about the local 
dynamical fluctuations is encoded in a purely local irre- 



ducible self-energy, S(w), entering the lattice Green func- 
tion for the model under consideration: 



G(k,uj) = G(e k ,uj) = 



1 



uj + p - e k - E(w) 



(5) 



To solve the model in d — oo requires a reliable way to 
solve the single impurity Anderson model (SIAM) em- 
bedded in a dynamical bath described by the hybridiza- 
tion fn. A (a;). There is an additional condition that 
completes the selfconsistency: 



/ 



deG(e,uj)p (e) 



1 



uj + fi — A(w) — S(w) 



(6) 



where po(e) is the free DOS (U = 0). The above eqns. (4)- 
(6) refer to the disorder-free Hubbard model. In d = oo, 
this is sufficient to compute the transport, because the 
vertex corrections in the Bethe Salpeter eqn. for the con- 



ductivity vanish identically in this limit 15 1. Thus, the 



conductivity is fully determined by the basic bubble dia- 
gram made up of fully interating local GFs of the lattice 
model. 

The optical conductivity is computable in terms of the 
full d — oo GFs as follows 15 1: 



a xx (iuj) = — I depo(e) V" G(e, iv)G(e, iv + iw) (7) 

IU1 



and the dc resistivity is just Pdc{T) — l/a xx (0,T). 

To solve the impurity model, we have used the iter- 
ated perturbation theory adapted to our model. We use 
the IPT because it is analytic in U and yields results in 
good agreement with exact diagonalization studies [^5| . 
In view of the ability of the IPT to reproduce all the 
qualitative aspects observed in <j xx {u>), we believe that is 
a good tool in the present case. We mention that we have 
extended the IPT off half-filling to finite temperatures (as 
far as we are aware, such an undertaking has not been 
carried out earlier). We have checked that our results 
agree fully with those obtained by others at T = 0. 

The photoemission lineshape is given simply by, 



iPES 



(uj) = p(ui - p,T)f(uj - p) 



(8) 



where f(ui) — [& ) ' T + 1) _1 is the Fermi fn. 

To make closer contact with experiment, we notice 
that, in CYO2, Gr exists in two crystallographically in- 
equivalent sites, resulting from a particular kind of or- 
bital ordering || . The potential fluctuation arising from 
this is of order of the crystal field splitting (~ 1 — 1.5 
eV), providing an intrinsic source of disorder scatter- 
ing (especially at T ~ 300K at which the PES experi- 
ment was performed) , further reducing the coherent spec- 
tral weight at p. Hence, we have considered the ef- 
fects of scattering off an additional local, static and ran- 
dom potential v = U/4, assuming a binary distribution, 
P(vi) = [S(vi) + 5(vi — v)}/2 (an additional assumption 
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is that the material is composed of an equal composition 
of the two inequivalent sites). To this purpose we have 
used the IPT (for correlations) in combination with the 
CPA for disorder, in a selfconsistent way | ft9f . 



0.60 




-4.0 -2.0 0.0 2.0 4.0 



(O/D 

FIG. 1. The density of states (DOS) for U/D = 3.0, 
S — 0.2, T = 0.017? and two values of local disorder po- 
tential: v = (continuous) and v = U/A (dotted line). The 
dot-dashed line is the DOS for the same parameters as above 
and v = (7/4 but with T = 0.173. 

In Fig. H we show the local spectral density, p(u>) 
for our model at different T for U/D — 3.0 This ratio 
was chosen because we are interested in dealing with the 
metallic phase here. It is clearly seen that the collective 
"Kondo" resonance is strongly T dependent because of 
particle-hole pair generation; above T ~ T co h = 0.13, 
the coherent feature is completely washed out, and the 
physics is understood in terms of carriers scattered in- 
elastically off effectively unquenched orbital moments, in 
analogy with the situation in the usual Hubbard model. 
This has important consequences for the optical and 
spectroscopic response, as we describe below. 

At low T, the optical conductivity shows a sharp Drudc 
peak characteristic of a correlated FL metal. This peak 
carries a small part of the total spectral weight, and 
most of it is distributed over a wide energy scale of 3eV, 
comprising the rather well pronounced mid-IR absorp- 
tion and higher energy features (Fig. ||). Increasing T 
(T = 0.173), we see that the Drude weight is reduced 
and the higher energy features are smoothened out, again 
resembling qualitatively the situation in CVO2 §• 

The d.c resistivity (not shown) is computed from the 
lu = part of the optical conductivity, and clearly ex- 
hibits the quadratic behavior in T: 

p dc (T) ~ T -\T) = 7r 2 ([/ 2 /73 3 )n/2(l - n/2)T 2 (9) 

With an effective mass enhancement of approximately 3 
for Uab/D = 3.0, n = 1.3.10 28 m~ 3 , p(ji) = O.tteV' 1 and 



n = 0.8, we obtain A/7 2 ~ 4 — 5 * 10 5 , as reported in 
Ref. §. 



0.60 
0.45 

3 

D S 0.30 
0.15 
0.00 









T=0.01D 






v=U/4 


T=0.01D 






v=U/4 


T=0.1D 




\ \ 








\ \ 
















\ 

\ V 








V \ 















0.0 1.0 2.0 3.0 

CO/73 



FIG. 2. Real part of the optical conductivity a xx for 
L7/73 = 3.0, 5 = 0.2 v = (continuous) and v = L7/4 (dot- 
ted), all at T = 0.0173. The dot-dashed line shows a xx for 
17/73 = 3.0, 5 = 0.2, v = 17/4 and T = 0.173. The computed 
result shows all the features observed in Ref. [9]. 

Next we turn to angle integrated photoemission, which 
measures the occupied part of the single-particle DOS, 
in Fig. ||. As expected from the T-variation of the DOS, 
we see a sharp Fermi step at T = 0.0173, a reduced step 
feature at T = 0.113, and completely incoherent response 
in the pseudogap phase. A comparison of the results with 
those observed experimentally shows that the calculated 
spectral weight at the Fermi surface (u> = 0) is quite 
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FIG. 3. The integrated photoemission lineshapes for our 
model with U/D — 3.0, 8 — 0.2 and for two values of the 
disorder strength v and the temperature T. Introduction of 
disorder brings the lineshape much closer to that observed ex- 
perimentally by moving low energy spectral weight to higher 
energy. 
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small, but still finite, while complete depletion of the 
low energy spectral weight is observed experimentally. 
The almost complete depletion of the low-energy spec- 
tral weight observed experimentally is then interpreted 
in terms of the T-dependence of the local spectral den- 
sity. If the temperature of measurement is higher than 
the lattice coherence scale, the FL-like feature is washed 
out, reducing the spectral weight at /z, and making the 
intensity resemble that of a semiconductor. In this con- 
text, we would like to mention that the PES experiment 
was carried out at T — 300 K, which is pretty high. The 
calculation carried out above does yield an incoherent re- 
sponse at T — 0.1D and the spectral weight at /i is quite 
small but is still non- vanishing at fi. 

As we have mentioned before, the results for a moder- 
ate disorder v — U/4 resemble experimental observations 
more closely (fig. 3). However, it may also be the case 
that short-ranged orbital correlations enhance the ten- 
dency to open up a pseudogap, reducing the DOS beyond 
that found in the d = co calculation. At this moment, 
this discrepancy between theory and experiment remains 
unclear. A more precise undertaking should also include 
the real bandstructure of CrOi', in this context, we no- 
tice that the LDA calculations already yield a dip in the 
DOS at Ep ||. We are of the opinion that the real situ- 
ation involves a combination of all the above factors, and 
plan to address these in a future publication. We suggest 
that more information concerning the T-dependence of 
the low-energy spectral weight, could be obtained from 
a simultaneous study of optical and PES spectra taken 
at different temperatures across the ferro-para transition. 
This could be used to discriminate between different sce- 
narios proposed for GYO2 ■ If the ideas proposed here are 
correct, one expects a T-dependent build-up of coherent 
spectral weight at low energy as T is lowered through the 
para-ferro transition. 

Finally, we observe that the model hamiltonian pro- 
posed in this paper can explain the orbital polarization 
observed by Stagarescu et al. ||. Indeed, for a generic 
choice of the t a ,a = 1,2, with t\ 7^ t2, an orbital po- 
larization is spontaneously generated in our model. A 
complete analysis of these effects requires consideration 
of actual bandstructure within the d = 00 approximation 
that we've used here, and we plan to treat this issue in 
the future. 

In conclusion, we have shown, based on the LDA+U 
results of Korotin et al. ||, how the various observed 
"strongly correlated" properties of the half-metallic fer- 
romagnet CrOi can be understood by invoking the role 
of dynamical orbital correlations and Hund's rule dou- 
ble exchange in the t2 g sector. The treatment presented 
here should also be valid for other transition metal based 
half-metallic ferromagnets H . Our results are consistent 
with the view that LSDA calculations || need to be sup- 
plemented by treatments including dynamical correlation 
effects to understand the physics of CrC>2 ■ 
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